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Abstract. The differential and total cross sections for the scattering of muonic, pionic, kaonic and antipro¬ 
tonic hydrogen in excited states from atomic hydrogen have been calculated for the purpose of atomic 
cascade calculations. The scattering problem is treated in a fully quantum mechanical framework which 
takes the energy shifts and, in the case of the hadronic atoms, the widths of the ns states into account. 
The validity of semiclassical approximations is critically examined. 

PACS. 34.50.-s Scattering of atoms and molecules - 36.10.-k Exotic atoms and molecules (containing 
mesons, muons, and other unusual particles) 


1 Introduction 


Exotic hydrogen-like atoms are formed in highly excited 
states, when negative particles ,Tr~, K~...) are stopped 
in hydrogen. The deexcitation of exotic atoms proceeds via 
many intermediate states until the ground state is reached 
or a nuclear reaction takes place. Despite a long history 
of theoretical and experimental studies (see and ref¬ 
erences therein) the kinetics of this atomic cascade is not 
yet fully understood. The current generation of experi¬ 
ments with exotic hydrogen-like atoms addresses a num¬ 
ber of fundamental problems using precision spectroscopy 
methods, the success of which relies crucially on a better 
knowledge of the atomic cascade. 

In the case of the laser spectroscopy of the Lamb shift 
in muonic hydrogen , the goal is to determine the proton 
charge radius with an accuracy of 10“^ from the energy 
splitting between the 2s and 2p states. This will remove 
the major theoretical obstacle in the precision calculations 
of the hydrogen Lamb shift, thus extending the limits of 
the most stringent test of QED in a bound system. The 
feasibility of this experiment depends on the population 
and the lifetime of the metastable 2s state of p~p, and a 
reliable model of the cascade kinetics is essential for this 
issue. The experiment on precision spectroscopy of pionic 
hydrogen is expected to determine the irN scattering 
length with a precision better than 1% by measuring the 
nuclear shifts and widths of the K X-ray lines. At this level 
of precision, the Doppler broadening corrections to the line 
width become important, and they must be reliably calcu¬ 
lated from a cascade model. In the precision spectroscopy 
of antiprotonic hydrogen 0, the Doppler broadening of 
the L X-ray lines must be taken into account when the 2p 
nuclear widths are determined from the X-ray line profile. 


The kinetics of atomic cascade is described by the mas¬ 
ter equation involving all significant processes with the ex¬ 
otic atoms (cascade mechanisms). The deexcitation mech¬ 
anisms include radiative. Auger, and Coulomb processes 
where the transition energy between states with different 
principle quantum number n is carried away mainly by 
photon, electron, and the recoiling particles (including the 
exotic atom itself), correspondingly. While the deexcita¬ 
tion processes are obviously essential for the atomic cas¬ 
cade, the role of the collisional processes preserving the 
principal quantum number n is not less important than 
that of the deexcitation. The Stark transitions nl nl' 
{V 7 ^ 1), affect the population of the nl sublevels. Together 
with the elastic scattering nl nl they decelerate the ex¬ 
otic atoms thus influencing their energy distribution dur¬ 
ing the cascade. In hadronic atoms, the role of the Stark 
mixing is especially important as it results in a strong ab¬ 
sorption during the cascade by feeding the ns states which 
have absorption widths much larger than the states with 
I > 0 . 

In the literature starting with the paper of Leon and 
Bethe , Stark mixing has often been treated in the semi¬ 
classical straight-line-trajectory approximation [|^, ^, |^ 
0 - Due to the broad use of this relatively simple model 
it is desirable to know its accuracy in comparison with 
more advanced and realistic quantum mechanical calcu¬ 
lations. By introducing phenomenological tuning param¬ 
eters in the Stark mixing rates one is able to repro¬ 
duce the measured X-ray yields and other experimental 
data. However, the ultimate goal of ab initio cascade cal¬ 
culations demands more accurate results for the collisional 
processes. 

A fully quantum mechanical treatment based on adia¬ 
batic potentials was given in [^2|-^. However, the shifts 
and widths of the ns states which become important in 
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the final part of the cascade were not included in this 
framework. The deceleration and radiative quenching of 
muonic hydrogen in the metastable 2s state were stud¬ 
ied in a close-coupling framework in Hi- We reexamined 
the same problem in p^ ] avoiding some of the approxi¬ 
mations used in . As the close-coupling model can be 
straightforwardly modified to include nuclear absorption 
in hadronic atoms ||^ , it is well suited for describing the 
collisional processes during the atomic cascade. 

In this paper, we present a unified treatment of Stark 
mixing, elastic scattering, and, in the case of hadronic 
atoms, nuclear absorption during collisions. For the time 
being, we restrict our calculations to exotic hydrogen¬ 
like atoms. The paper is organized as follows. In sect. || 
we present the quantum mechanical close-coupling frame¬ 
work used for the calculation of the scattering of x~p 
atoms in excited states from atomic hydrogen. The same 
processes are treated in the semiclassical approximation 
in sect. The results (differential, partial wave and total 
cross sections) are discussed in sect. ^ and summarized in 
sect. ||. 

Unless otherwise stated, atomic units {Ti = e = mg = 
1) are used throughout this paper. The unit of cross sec¬ 
tion is ttg = 2.8 • cm^, where og = ^ ^‘2 is the 

electron Bohr radius. 


2 Close-coupling calculation of the cross 
sections 

In this section, the close-coupling framework developed in 
ref. 1^ for the 3-body reaction 

{X~p)nl -b H ^ {X~p)nl’ + H (1) 

will be generalized to include absorption effects in hadronic 
atoms. The following notations are used: the negative par¬ 
ticle x~ with mass rux and the proton with mass mp form 
an exotic atom with the total mass M^p = rrix + rup 
and the reduced mass = mx'mp/Mxp. The coordinates 
used in the calculations are explained in fig. |^. The rela¬ 
tive orbital angular momentum of x~p and H is denoted 
by L, the internal x~p orbital angular momentum by 1, 
and the total orbital angular momentum by J = L -|- 1. 


2.1 The effective x p — H interaction 


The exotic atom x p is described by the Hamiltonian 


^xp - 



- + AV 

r 


( 2 ) 


where AV includes all effects beyond the standard non- 
relativistic Coulomb two-body problem: vacuum polar¬ 
ization, finite size effects, and, in case of hadronic atoms, 
strong interaction between the two particles. For the con¬ 
sidered systems, AV can be treated as a perturbation re¬ 
sulting in the shift AEni of the Coulomb energy levels: 


{nlm\Hxp\nlm) = E^i 




+ AE, 


'nl 


H 



Fig. 1. Coordinates used for the effective 3-body system x p— 
H: R is the vector from the target proton to the center of 
mass of the exotic atom, r is the relative vector of the x~p 
system. The solid angles of R and r are denoted by U and ui, 
respectively. 


where the ket \nlm) denotes the nlm state of the stan¬ 
dard Coulomb problem. In the case of hadronic atoms, 
the atomic cascade is often terminated before reaching 
the ground state due to nuclear reactions like 

TT~p — >■ 7r°n, yn , (4) 

K-p^ E^TT^,E°TT^,An . (5) 

This effect is described by the imaginary part of the com¬ 
plex energy shift AEni 


lm{AEni) = -r„,/2 (6) 

where the width Eni is the nuclear reaction rate from 
the nl state. For the collisional processes with the exotic 
hydrogen-like atoms, it is sufficient to take into account 
only the widths of the ns states. The n dependence of the 
hadronic part of the complex nuclear shift is described by 
the Deser formula ||^ 


A ZT’had 

;^had _ ^-^Is 


AE'^T = 


( 7 ) 


The Hamiltonian for the x p — H system is given by 




H. 


xp 


( 8 ) 


where 


= 


MxpM]i 


M, 


xp 


■ ATh 


( 9 ) 


is the reduced mass of the [x~p + W) system, and U(r, R) 
is the Coulomb interaction of the exotic atom with the 
hydrogen electric field. In most cases of cascade studies, 
the effects of hydrogen excitation or ionization in the col¬ 
lisions without changing n can be neglected. The electric 
field of a hydrogen atom in the ground state has the form 


F(R) = |f(R) 


( 10 ) 


where 


E{R) = —(1 + 2RV 2i?2)e-2« 


( 3 ) 


( 11 ) 
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The potential energy is a function of only three variables 

V{v,^)=V{z',t',R) (12) 

where r' = {x\y\z') is r in a rotated coordinate system 
with the z'-axis taken along R. The electric field of a hy¬ 
drogen atom is sufficiently strong to mix the I sublevels 
of the x~p atom during a collision at distances of a few 
atomic units, oq, which are much larger than the size of 
the x~p at low n. Since Stark mixing is essentially a long¬ 
distance process, one can use the dipole approximation for 
the potential 

V{r,R) = z'F{R) . (13) 

This can be easily generalized to scattering from other 
target atoms: one uses the electric field generated by the 
nucleus and the electronic charge density. 

Given the above defined interactions, we solve the time 
independent Schrodinger equation 

R'0(r, R) = if'^(r,R) (14) 


The ket |n/yl)' denotes the eigenstates of the Coulomb 
problem in the rotated coordinate system. 

Because of rotational invariance the quantum num¬ 
bers J and M are conserved and the radial wavefunc- 
tions (jijiA^R) are independent of M. The expansion ( 0 ) 
is convenient for computing matrix elements of the po¬ 
tential V(z', r', R). In the dipole approximation, the non¬ 
vanishing matrix elements of z' correspond to |Z\/| = 1 
and AA = 0 where one has 


\nlA\z'\n{l — l)A)' = 


2>n j {P — A^){rF — P) 

(2; + i)(2;-i) 


The basis states (|l^ are not the eigenstates of but 
the matrix elements of \P can be easily obtained by using 


= (J - 1)2 = J2 + i2 _ l'_J+ . (20) 

Together with the results and notations of appendix ^ 
this gives (see also ref. [pOf) 


with the standard boundary conditions for this multichan¬ 
nel scattering problem using the close-coupling approxi¬ 
mation. The internal x~p wave function is expanded into 
the set of rP Coulomb wave functions with the same quan¬ 
tum number n: 

V'(r,R)=R-i ^ </>jM(i?)^i*M(^?)x«M(r') . (15) 

JMIA 

The functions XniAiP) are the Coulomb wave functions 
in a rotating coordinate system with the z'-axis chosen 
as the quantization axis and are the correspond¬ 

ing rotation functions (see appendix ^). The functions 
XniA{P) are related to the space fixed functions Xn/m(r) = 

rni{r)Yi 

m (w) by 


I 

Xnim(r) = Dl^^{n)xuiA{r') . (16) 

A=-l 


The expansion (^) leads to the set of coupled sec¬ 
ond order differential equations for the radial functions 
4>JIa{R) 


1 d^ 


+ Enl — E 


^(I>jia(R) 


E 


2pL dR2 

n]JMA'V\lP\n]JMAl) 


2/xR2 


\nl'A'\V{z\r',R)\nlA)']4>ji,A'{R)=0 . (17) 


The basis states \n; JMAl) are simultaneous eigenstates 
of El^p, J2, and P with eigenvalues E^i, J{J +1), 

M, A, and l(l + 1), respectively, and are given by 


(n; JMA'l'l L^ln;JMAl} = 

SwSaa' (j(J+1)+ ?(; + !)-2712) 

— Sw ^l5yl+lyl/A_(J, A)X+{1, A) 

+i^yi-iA'A+(J, 7l)A_(/, yl)) . (21) 

The rotated basis functions were used by Carboni and 
Fiorentini to study (/r“p) 2 s+H collisions in an approx¬ 
imation where the terms 6wSaa>{1{1 + 1) — 271) in eq. d^ ) 
were neglected. 

One can get a partial decoupling of the equations 0 ) 
by using the following expansion 

V^(r,R)=i?-i ^ (ju{R)yi^in,co)r^i{r) (22) 

JMLl 

where the functions 

= Y {LlMLm\JM)YLMd^)Yi^{io) (23) 

MLm 


are simultaneous eigenfunctions of J2, L2, C, and Jz with 
eigenvalues J(J-l-l), L(L-l-l), l{l + l), and M respectively. 
The system of the radial Schrodinger equations for the 
functions ^juiR) has the form 


7 1 , LjL + l) 

V 2 fi di ?2 2 /rR 2 

+ Y{n;L'l'JM\Viz',E 


+ Enl — E'j^jLi{R) 


,R)\n;LlJM)^jLn'{R)=0 


L'l' 


(24) 


with the basis states 


|n; JMTlO = |JM7l)|n;7l)' , 

{Q\JMA) = y^^7?](;n(^2) 
(r'|n/7l)'= x„iyi(r') • 


(18) {n,v\n-LlJM)=yif^{n,uo)rni{r) . (25) 

Due to the parity conservation, the value P = (—1)^+* 
is conserved, and, as a result, the differential equa¬ 
tions ( 0) are decoupled into two sets of n{n + l)/2 and 
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n(n — l)/2 coupled equations for P = 1 and P = —1 cor¬ 
respondingly. The systems of the equations ( 0 ) and 
are related to each other by the linear transformation 


1 TT 


(t>JiA{R) = ^ y^^UALpiLiR) 


AA' 


2.2 Cross sections 


(27) 


dcr„i 


>7ll' 


1 


k' 


df? 


{21 + 1) k 


El/"' 


m^nV m' 


(30) 




{21 + 1 ) fc 2 
X E \{n\L'l'JM\S-l\n]LlJM)\^ 


(26) 


where the coefficients are given by eq. ( |A10| ). The ma¬ 
trix elements of the potential energy can then be obtained 
from those of the rotated basis by using the coefficients 
^AL 

{n-,L'l'JM\ V{z',r',R)\n]LlJM) = 

'^u'a'L'Ual \nl'A'\V{z\r',R)\nlAy . 


JMLL' 

1 TT 

( 2 / + 1 )F 


E(( 2 J+ 1 ) 


X -Mri;LlJM)f'^ .(31) 

LL' 

The corresponding transport cross sections are given by 

, d(Jnl—*nl' 


= / dl 2 (l-cos 6 »)- 


df} 


(32) 


The scattering matrix was calculated numerically using a 
version of the variable phase method described in ap¬ 
pendix In order to treat the ns states of hadronic atoms 
as normal asymptotic states the absorptive term (^) was 
switched off for the distances between x~p and H larger 
than Soq. The absorption from the ns states between the 
collisions can be easily taken into account by means of a 
cascade model. 

The use of the dipole approximation in the quantum 
mechanical framework makes it necessary to introduce the 
regularization parameter Pmin as explained in appendix ^ 
The dependence of calculated cross sections on Pmin will 
show how sensitive the results are to the short distance 
behavior (we will show a few examples in sect. 

The scattering amplitude for the transition nlm 
nVm' is given by 

/ nlm^nl'm'{^) = /— E Yl'Mi{^) 

2zVfc k ^ 

X (n; L'l'M^rn'lS - l|n; LZ0m)y£o(0,0)) (28) 

where f2 is the CMS scattering angle, k and k' are the 
CMS relative momenta of the initial and final state corre¬ 
spondingly. The S-matrix elements in ( p^ are related to 
the matrix elements between the basis states ( [ 1 ^ ) by the 
relation 

{n] L'I'\S\n] LlQm) = 

{L'I'M^m!\Jm){Jni\LlQm){n\ L'l'Jni\S\n] LlJm) 

j 

(29) 

The differential and total cross sections for the transitions 
nl nV are given by 


In the case of hadronic atoms, the scattering matrix is not 
unitary because of the absorption. The cross sections for 
the absorption processes are given by 

Cni^abs = p E 
j 

Xll- ^ 


(1-72^ E \{n-,L'l'JM\S\n-LlJM)y)) 

^ ’ LL'l' 


(33) 


The differential cross sections (po| ) and the absorption 
cross sections (^lare used in the detailed cascade models 
as described in . When less detailed information is suf¬ 
ficient, Z-average cross sections defined below can be used. 
In particular, in those cases where the rates for collisions 
without change in n are much larger than other cascade 
rates, the approximation of the statistically weighted dif¬ 
ferential cross section is useful: 


du n—a 

Afl 


= iE( 2 ^ 


n 


1 ) 


du. 


nl- 


di? 


(34) 


As a measure of the overall strength of Stark mixing, one 
can use the Z-average Stark cross section: 


cst = 


4 e 

i^v 


{21 -\- 


(35) 


For hadronic atoms with the strong absorption in the 
s states, we define the statistically weighted differential 
cross section (34) and the average Stark cross section (p|) 
to include only terms with I > 0 and I' > 0. The average 
absorption cross section is defined by 


1 


<^abs — 


n^-1 


E^( 2 Z + l)o'nl- 


>-abs 


(36) 


1^0 


which gives a measure of the absorption strength under 
the assumption of statistical population of the I 7 ^ 0 sub- 
levels. Strong absorption can also take place between the 
collisions (if the hadronic atom leaves the collision zone in 
an s state) which is not reflected by eq. (p^. To this end 
we define the average cross section 


n2 - 1 

1^0 


{21 - 1 - l)ani- 


2 


(37) 
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which describes the transitions to the ns state from the 
statistically populated nl sublevels with I > 0. Whether 
the s state is completely or partially depleted between 
the collisions depends on the type of atom, the quantum 
number n, the density of the target, and the kinetic energy. 
To estimate the upper limit of the absorption, we define 
the maximum absorption cross section as the ^-averaged 
sum of the cross sections for nuclear absorption during 
and after collision 


^max abs — *^abs Ua 


(38) 


The Stark mixing, deceleration, and absorption rates, 
which are often used in cascade calculations, are defined 
by the formulas 


Ast = Nvast , 

, _ „ MiiMxp 

'^dec — 


(Mh + M,p)2 
Aabs — ^ abs 


Nvut, 


(39) 

(40) 

(41) 


where N is the target density and v the velocity of the 
exotic atom. When a significant part of the nuclear reac¬ 
tions takes place between the collisions, the absorption is 
better described by the effective absorption rate defined 
as following 


Aeff abs — Aabs 4“ 


'av —>-ns 


1 + ^ 


■*nl 




(42) 


semiclassical treatment applicable to many cases of prac¬ 
tical interest. The model that has been used most often 
is the straight-line-trajectory approximation Sim 
where the small neutral exotic atom is considered as mov¬ 
ing along a straight line with constant velocity through 
the electric field of the target atom]^. The time depen¬ 
dent electric field causes transitions among the sublevels 
of the exotic atom, which are treated quantum mechani¬ 
cally. This approach was usually used for the calculation 
of the Stark mixing rates, but, as discussed below, differ¬ 
ential and absorption cross sections can be calculated as 
well. 

A semiclassical description of our scattering problem 
is obtained by treating some of the 6 variables (R and 
r) as classical time dependent variables. The remaining 
variables correspond to the quantum mechanical part of 
the system that is described by the wave function ipit) 
satisfying the Schrodinger equation 

(44) 

where H®®(t) depends on t through the classical variables. 
The wave function ijj(t) is expanded into a set of orthonor¬ 
mal basis states 

fj{t) (45) 

3 

leading to the time dependent Schrodinger equation 


In the case of very strong absorption during the collisions, iaj{f) — ^ H^^{t)ak{t) (46) 

the relation ‘C Aabs holds, and therefore j. 


Aeff 


abs 


'^abs (for 




(43) 


When absorption from the p states is important, the effec¬ 
tive rate for p state absorption is defined analogously to 
eq.(42) by considering statistically populated I > 1 states. 
A simple comparison of the energy dependent rates for 
the different processes cannot substitute detailed cascade 
calculations using the detailed cross sections eqs. (^,|3|) 
but may be helpful for getting a quick overview. We shall 
present a few examples in sect. 0. 


3 Semiclassical approximation 

As the number of coupled second order differential equa¬ 
tions in the quantum mechanical model of sect. ^ grows 
as n^, a simpler framework is desirable for high n states. 
If the collision energy is sufficiently large, one can expect 
that the relative x~p — H motion can be treated classi¬ 
cally. A rough estimate of the minimum kinetic energy T 
for which a classical-trajectory description is valid can be 
obtained for the requirement that a large number of par¬ 
tial waves L ^ 2aok (2ao being the approximate range of 
the interaction) contribute to the cross section; that gives 
for the kinetic energy of muonic hydrogen T > 0.7 eV 
at L > 10. As known from experiment (see Q] and ref¬ 
erences therein), the exotic atoms can reach kinetic ener¬ 
gies of several eV during the cascade, and this makes a 


which must be solved with appropriate boundary condi¬ 
tions. 

To establish a connection to results in the literature we 
will first discuss the simple fixed field model of Leon and 
Bethe 00 where the R motion is assumed to be clas¬ 
sical and the x~p atom is treated quantum mechanically. 
The assumptions are as follows: the x~p moves along a 
straight line with constant velocity v {R{t) = ^{vtY + 
where p is the impact parameter), only transitions within 
the states with the given principal quantum number n 
are considered, all the states are degenerate, and the 
electric field from the target atom is directed along the 
quantization axis of the x~p. After expanding the inter¬ 
nal wave function of the exotic atom into the Stark eigen¬ 
states (|nniA), m = 0, ...,n — |A| — I) one is left with a 
single channel scattering problem 

ib^nniA(f) — (7?(t))u^yj,^yl (t) (47) 

where 

S?7 

VnnMR) = i^ - i2n,-n+\A\ + l)FiR) . (48) 

Equation ( ^^ is solved with the boundary condition 

anniA(-Oo) = 1 (49) 

^ Another possibility, which we will not consider here, is to 
use deflected trajectories. 
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for a range of values of the impact parameter p 


where is a x matrix given by 


= exp J V{R{t))dt^ 


(50) = 


sc/,\ L{L + 1) — J{J + 1) 


2pR^{t) 


+ V{R{t)) + AE . (58) 


The eikonal phase shift function 

nOO 

Xip) = - / V{R{t))dt (51) 

J oo 

is used (we take J+1/2 = kp, where k = p,v is the relative 
momentum) to obtain the scattering amplitude 


^ E (2J + - l)Pj(cos0) 


j=o 


and the differential cross sections 


J —eikonal 
^^nniA 

df2 


= I f 

\J 1 


eikonal 

nniA 




(52) 


(53) 


This m odel is not sufficiently accurate for our purposes. 
In sect. ^ we generalize it to include different thresholds, 
nuclear absorption during collisions, and correct angular 
coupling between the substates. 


3.1 Semiclassical model 


The potential matrix has the elements 

V,j{R) = {n]L'l'JM\V\n-LlJM) (59) 

as in eq.(^J). 

The eqs. ( |5^ ) are integrated from —t„iax to tmaxj where 
tmax is chosen so large that the potential can be neglected 
for distances larger than i?(tmax), with the boundary con¬ 
ditions 

^(-tmax) = I ■ (60) 

The semiclassical scattering matrix is given by 

S'®® = QA(tmi,^)Q (61) 


where Q is the diagonal matrix given by 

. L(L -l-l) — T(T + 1) 


Q = exp (^i 


2pk 

-|-zRe(Z\ii/)tniax^ 


arctan(vti„ax/p) 


(62) 


With the semiclassical scattering matrix and the for¬ 
mulas from sect. results for differential, total and ab¬ 
sorption cross sections can be obtained. 


In this case only the radial R motion is considered to be 
classical, the other five variables (17, r) are kept quantized. 
The problem of x~p scattering from hydrogen becomes 
a multichannel scattering problem with different channel 
momenta and orbital angular momenta. We collect the 
(complex) energy shifts and the angular part of the kinetic 
energy of the different channels in the diagonal x 
matrices AE and L, respectively. The collision is specified 
by the CMS collision energy 

Ecm = ^ (54) 


3.2 Fixed field model 

The semiclassical approximation simplifies the numerical 
calculations: instead of a system of second order differen¬ 
tial equations the same number of first order differential 
equations must be solved. A further simplification can be 
made by neglecting the coupling between the internal x~p 
angular momentum 1 and the orbital angular momentum 
L. This corresponds to the approximation 

L = J (63) 


and the angular momentum J. Neglecting the deflection 
of the neutral x~p atom we take the classical motion to 
be 

R = R{t) = ^(ut)2 -k p2 (55) 

where p = J{J + ^)/k. This introduces some ambiguity 
into the model because it requires a common motion R{t) 
for all channels. In this paper we take the common mo¬ 
mentum k to be that of the I = {n—1) states for hadronic 
atoms while we use the ns state in the case of muonic 
hydrogen. 

We expand the quantum mechanical part of the system 
into the states 

\aj) = \n;LlJM) (56) 

and obtain the time dependent Schrodinger equation in 
matrix form 

(57) 


in the case discussed above in sect. 3.1. In this approxi¬ 
mation, the quantum number A is conserved in addition 
to J and M (see fig. ||). 

By expanding the solution into the rotated basis states 
{\aj) = |n; JMAl)) one finds the time dependent Schrodinger 
equation in matrix form 


iA(t) = {ZF{R(t)) + AE) A(t) , (64) 

A{-tmax) = I (65) 

where 

Zij = '{nlA\z'\nl' A)' . (66) 

Equation (|^) must be solved for each value of A, |A| < 
n. If the states with I > 0 are taken to be degenerate, 
eq. (|^ with A ^ 0 decouple completely in parabolic co¬ 
ordinates 1^,^ and is easily integrated as shown above. 
This is also the case for A = 0 when the ns energy shift 


iA{t) = H^^{t)A{t) 
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(a) 


L-J 


(b) 

A 


< 



2 

1 

0 

-1 

-2 


2 

1 

0 

-1 

-2 


1 = 0 1 2 


1 = 0 1 2 


Fig. 2. Coupling between the sublevels in an —H collision. The correct coupling (a) conserves parity. In the fixed field model 
(b) transitions are only possible between states with the same eigenvalue A of x~p angular momentum along the interatomic 
axis. 


is negligible. If not, one must solve n coupled first order 
equations. 

The fixed field scattering matrix is defined by 

(67) 

where the diagonal matrix QFF jg given by 

QFF = exp (^iRe(AE)tmax'j ■ (68) 


The fixed field scattering amplitude is given by 

fnlA — ynl' a{^') — ^ 

X ^(2J + I)(n; JMAl'\S^^ - l|n; JMAl)Pj{cos 9) 

J 

(69) 


and the differential cross section by 


dI2 


2Z + 1 


^ ^ \ fnlA- 


'll' A 


(0)1^ 


(70) 


The matrix elements in the r.h.s of eq.(|6^ are actually in¬ 
dependent of M since this quantum number is conserved. 
The cross sections for the processes nl nV and nl 
absorption are given by 

I TT 

(Tnl^uV = 

X^(2J+1)^ \{n-JMAl'\S^^-l\n]JMAl)\^ 

J .4 

(71) 

and 

CTni^abs = 2/ I + 1 )(( 2 ^ + 1 ) 

-^\{n-JMAl'\S^^\n;JMAl)\^^ . (72) 

Al' 


4 Results 

Using the methods described in sects. and ^ we have 
calculated the cross sections for the collisions of the 
TT~p, K~p, and pp atoms in excited states with hydro¬ 
gen atoms. Our calculations had two major goals: first, to 
provide comprehensive sets of the collisional cross sections, 
which are necessary for detailed cascade calculations, and, 
second, to investigate the range of validity of the approxi¬ 
mate methods based on the semiclassical model and often 
used in the literature. The numerical calculations have 
been done for the principal quantum numbers n and the 
atomic kinetic energies that are of interest for the cas¬ 
cade calculations. The quantum mechanical framework of 
sect. H was used for the lower excited states n = 2 — 5, 
and the semiclassical calculations were done for the range 
of n up to n ~ 10. 

As the number of the calculated differential cross sec¬ 
tions is quite large (about 1200 for pL~p) only a small part 
of the results can be shown here, as we describe some 
main features of the calculated cross sections illustrating 
them with particular examples for different exotic atoms. 
Concerning the detailed results, they have all been used 
as input for the Monte Carlo kinetics code [^,|2j], and 
the results of the cascade calculations will be published 
elsewhere. 


4.1 Muonic hydrogen 

The muonic hydrogen scattering is the least complicated 
case because there is no nuclear absorption in the inter¬ 
action. The differential cross sections are known to have 
a characteristic shape with a strong forward peak and a 
pattern of maxima and minima as expected for 

the interaction that is essentially of a dipole-like type. 
Figure shows an example of the differential cross sec¬ 
tions for the elastic scattering 5s ^ 5s and Stark transi¬ 
tions 5s —> hp and 5s ^ 5g. The elastic cross section has 
a strong peak at zero scattering angle, while the Stark 
transitions reach their maxima at Hnite scattering angle. 
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Fig. 3. The differential cross sections for (/i~p) 5 s + H —> 
(l-i~p)sa,p,g + H vs. CMS scattering angle 6 at the laboratory 
kinetic energy T = 3 eV. 



Fig. 4. The statistically weighted differential cross sections 
for {p~p)n =5 + H —> {iJ,~p)n .=5 + H VS. CMS scattering angle 6 
at the laboratory kinetic energy T = 3 eV. The fully quantum 
mechanical (QM) results computed with JJmin^O.OS, 0.1, 0.15 
and 0.20 are shown with solid and short-dashed lines, the result 
of the semiclassical (SC) model is shown with a dashed line, 
and that of the hxed held model (FF) with a dash-dotted line. 

The peaking in the forward hemisphere is much less pro¬ 
nounced for larger changes in quantum number 1. 

To compare the results of different methods from the 
viewpoint of practical applications it is better to look at 
cross sections averaged over some appropriate distribution 
over I or kinetic energy T since many tiny details will be 
washed out anyway in the cascade evolution. For the pur¬ 
pose of illustration as well as for simple estimates, the sta¬ 
tistically weighted cross sections are especially useful. Fig¬ 
ure H shows an example of the statistically weighted differ¬ 
ential cross section for n = 5 calculated in the fully quan¬ 
tum mechanical model for four values of the cut-off pa¬ 


rameter Rmin (0.05, 0.1, 0.15 and 0.2) in comparison with 
the semiclassical and fixed field approximation. While the 
cross sections for the individual transitions nl —>■ nl' are 
somewhat sensitive to the short range cut off the 

cut-off dependence smoothes out in the averaged cross 
sections. The semiclassical results are in a good agreement 
with the quantum mechanical ones for energies above 1 eV 
and n > 2. The fixed held model provides, on average, a 
fair agreement with the more accurate methods for the 
scattering. 

Another illuminating way to inspect the complicated 
structure of the differential cross sections is presented in 
fig- I that shows the partial wave cross sections for the 
reaction {fJ,~p) 5 s -b H —> {p~p) 5 p -b H at T = 3 eV. The 
partial waves can be divided in two groups corresponding 
to the regimes of ”weak-coupling” or ’’strong-coupling” 
behavior. The higher partial waves can be reliably de¬ 
scribed (J > fca « 18 for this example) in the semiclassi¬ 
cal approximation, with the partial cross section showing 
a smooth dependence on the total angular momentum J. 

For the lower partial waves, all I states are strongly 
mixed with each other, and the partial cross sections dis¬ 
play a strong dependence on both J and T. While the 
semiclassical approximation is not applicable in this sit¬ 
uation for individual amplitudes, it still makes a reason¬ 
able estimate for the average partial cross sections as they 
are mainly determined by the statistical weight of the fi¬ 
nal states. As long as the largest contribution to the to¬ 
tal cross section comes from the total angular momenta 
corresponding to the semiclassical regime, the semiclassi¬ 
cal approximation is adequate for all practical purposes. 
These two regions of J are also different with respect to 
the dependence of the Stark cross sections on the change of 
the orbital quantum number 1. In the semiclassical regime, 
the transition amplitude rapidly decreases with increasing 
change in (. In particular, the corresponding partial cross 
sections for the transition 5s ^ 5p are larger than the 
ones for the 5s —*■ 5d transition. 

Figure || shows the statistical average partial cross sec¬ 
tions for Stark mixing in the states n = 3,4, 5 in compari¬ 
son with the unitarity limit. In the strong coupling regime, 
a simple estimate for the average cross section can be ob¬ 
tained by assuming that the scattering phases are rapidly 
changing, so that they appear as being ’’random” (in the 
old picture of the Stark phase accumulated along a tra¬ 
jectory with a small impact parameter it corresponds to 
the so-called complete mixing when the initial state is for¬ 
gotten after the collision). For the higher partial waves, 
the average cross sections are limited not by the unitarity 
constraint, but by the centrifugal barrier which becomes 
so strong that it prevents the exotic atom from getting 
close to the hydrogen atom. 

The dependence of the total cross sections on the ki¬ 
netic energy is shown in fig. ^ for the transitions with dif¬ 
ferent change in the fj,~p orbital quantum number If — U = 
0,1,4: {p~p)^s -b H ^ {p~p)bs,p,g + H. As it was said be¬ 
fore, the semiclassical calculations are in a good agree¬ 
ment with the quantum mechanical ones for kinetic en¬ 
ergies above I eV. Below 1 eV, where only a few partial 
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Fig. 5. Partial wave cross sections for + H ^ 

(/i“p) 5 p + H vs. total angnlar momentnm J at the laboratory 
kinetic energy T = 3 eV. The fully quantum mechanical re¬ 
sults are shown with bars, the dashed and dash-dotted lines 
correspond to the semiclassical and the fixed field models. 



Fig. 6. Total angular momentum dependence of the average 
Stark cross sections for n = 3,4, 5 at the laboratory kinetic 
energy is T = 3 eV. The light solid line is the unitarity limit, 
and the light dashed line is the result for the statistical mixing. 


waves contribute, the agreement is still fair after averag¬ 
ing over some energy range. The fixed field model is in a 
fair agreement with the other two. The cross sections of 
the fixed field model tend to oscillate more than the cross 
sections computed in the other two models. For example, 
in the fixed field model the 5s state is coupled only to one 
of the 5g substates and the corresponding transition is 
described by one phase shift for each partial wave J. The 
quantum mechanical and the semiclassical models connect 
the 5s state with 5 of the 9 substates 5g (see fig. ||) for a 
given angular momentum J, and the average cross section 



Fig. 7. Cross sections for (^ p) 5 s T H —^ H vs. 

laboratory kinetic energy T. The fully quantum mechanical 
results (solid lines) are shown in comparison with the results 
of the semiclassical (dashed lines) and the the fixed field model 
(dash-dotted lines). 



Fig. 8. The energy dependence of the /-average Stark (solid 
line) and deceleration (dashed line) rates for (p“p )„^5 at 
15 bar in comparison with the radiative 5p —» Is rate (dot¬ 
ted line), the rates for n — 5 —> A external Auger effect ||^ 
(light short-dashed line) and the n = 5 ^ 4 Coulomb deexci¬ 
tation (light dash-dotted line). 


is smoother as it is distributed over a larger number of 
individual contributions. 

A brief overview of the Stark mixing and the deceler¬ 
ation in competition with the deexcitation mechanisms is 
presented in fig. || for a typical example of the g~p state 
n = 5 in hydrogen gas at 15 bar corresponding to 0.018 of 
liquid hydrogen density. The muonic hydrogen atoms ar¬ 
riving at the n = 5 state during the atomic cascade will, 
on average, undergo a few transitions changing the orbital 
quantum number and loose about half of their kinetic en- 
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ergy before Auger deexcitation or, less probably, the ra¬ 
diative transition takes place (assuming no other effects, 
like muonic molecule formation). Because the evolution of 
the kinetic energy distribution is important, detailed ki¬ 
netics cascade calculations are needed to treat this 

problem. 


4.2 Pionic hydrogen 


The scattering problem for hadronic atoms is more com¬ 
plicated than for the muonic atoms because of the nu¬ 
clear absorption and larger energy shifts. The previous 
fully quantum mechanical calculations based on adiabatic 
potentials ||^, ^ did not take the shifts and widths of 
the s states into account, and all the other studies were 
based either on the approximation suggested in Q or on 
the time dependent Stark mixing along classical trajecto¬ 
ries ^ . In this section we present the first results 
of the quantum mechanical calculations in comparison 
with the traditional approximations. The strong interac¬ 
tion shift and width of the Is state of pionic hydrogen from 
the final analysis of the PSI experiment were used 

to calculate the complex energy shifts of the ns states: 


Zll^ns — -^- H e, 


vp 

ns 


e'lf =-7,11 eV 


(73) 

(74) 


Fu = 0.87 eF 


(75) 


wheree^P is the energy shift due to the vacuum polariza¬ 
tion ||] 


= -3.24 eV , (76) 

e^P = -0.37eV , (77) 

e^P = -0.11eV . (78) 

The collisions with transitions between the states I > 0 
{n > 2) are qualitatively similar to muonic hydrogen with 
respect to Stark mixing and differential cross sections. The 
transitions to the ns states are less probable due to the 
strong interaction energy shift and the nuclear reactions 
taking place during collisions. 

To illustrate the effect of the complex energy shift we 
consider the cross sections a 2 p^ 2 p and the maximum ab¬ 
sorption cross section CTmax abs = cr 2 p^abs + Cr 2 p^ 2 s at 
r = 3 eV for different unphysical values of the 2p — 2s 
energy difference and the widths using the models de¬ 
scribed in sects. || and |[ Figure j^shows the dependence 
on the energy shift for the physical value of the width 
/ 2 s =0.11 eV. All three models feature a strong influence 
of the energy splitting |/? 2 p — Re(if 2 s)| on the Stark mix¬ 
ing with the 2s state: the nuclear absorption is much more 
likely when the energy splitting is small in comparison 
with the characteristic Stark splitting in the electric field 
of the target atom. For small \E 2 p — Re(if 2 s)|, the semi- 
classical model agrees well with the quantum mechanical 
results for both cross sections, but starts to deviate when 
the energy splitting is increased. 



E,p-Re(E, 3 ) (eV) 


Fig. 9. The pionic hydrogen cross sections for the elastic 
scattering 2p —> 2p and the maximum absorption from the 
2p state vs. the 2p — 2s energy difference at the laboratory 
kinetic energy T — 3 eV. The fully quantum mechanical results 
are shown with solid lines. The semiclassical results are shown 
with dashed lines and those of the fixed field model with dash- 
dotted lines. The light vertical line shows the physical value 
for the energy splitting. The 2s width is Tbs = 0.11 eV. 


The dependence of the cross sections on the T 2 s with 
the energy difference fixed at the physical value 1.26 eV 
(the sum of the strong interaction shift and the vacuum 
polarization) is shown fig. |^. For T 2 s > 0.5 eV, the 
absorption cross sections calculated in the semiclassical 
model are in good agreement with the quantum mechani¬ 
cal result, while the result of the fixed field model is about 
20 % lower. The good agreement between the semiclassi¬ 
cal and the quantum mechanical model can be explained 
as follows: when the nuclear reaction rate is high the ab¬ 
sorption process takes place immediately from a mixed 
2s — 2p state and is not very sensitive to the kinemat¬ 
ics of the 2s channel (the available phase space), which 
is treated incorrectly in the semiclassical approximation. 
Both the semiclassical and the fixed field models break 
down for smaller T 2 s, including the physical value, where 
an accurate treatment of the kinematics is necessary. 

In general, the fixed field model underestimates the nu¬ 
clear absorption during collision in cormarison with the 
semiclassical model. As shown in fig. g, the fixed field 
model allows only one of the 2p substates to be mixed 
with the 2s state and undergo nuclear absorption, whereas 
the correct parity conserving angular coupling used in our 
semiclassical model mixes two of the 2p substates with the 
2s state (in the coupled basis). For example, for the phys¬ 
ical values (/? 2 p — Re{E 2 s) = 1-26 eV and T 2 s = 0.11 eV) 
the maximum absorption cross section at 3 eV is increased 
by about 18 % when the semiclassical model is used in¬ 
stead of the fixed field model. 

In our quantum mechanical model, the interaction be¬ 
tween the exotic atom and the target hydrogen atom is ap¬ 
proximated by the dipole term, which is adequate for large 
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Fig. 10. The pionic hydrogen cross sections for the elastic 
scattering 2p 2p and maximnm absorption from the 2p state 
vs. the nuclear width F 2 s at the laboratory kinetic energy T = 
3 eV. The fully quantum mechanical results are shown with 
solid lines. The semiclassical results are shown with dashed 
lines and those of the hxed held model with dash-dotted lines. 
The light vertical line shows the physical value for the 2s width. 
The 2p — 2s energy splitting is 1.26 eV. 



R(a„) 


Fig. 11. The effective potentials for the system (tt p)n =3 + H 
with total angular momentum J = 4. The curves are labeled 
with the quantum numbers (ni,yl). 


distances. When the distance becomes small, this approx¬ 
imation breaks down together with the other ones (the 
close-coupling expansion into the basis of atomic states 
with the same n without taking into account symmetry 
requirements for identical particles). Neglecting rotational 
coupling and energy shifts, the effective potential energy 
in the dipole approximation for 7r“p — H system is given 

by 

Kff.(-R) = -z - {2ni-n+\A\ + l)F{R)+ ^ ^ (79) 

2pR^ 


Fig. 12. /-average Stark, transport, and maximum absorption 
cross sections for (7r“p)„=3 -|- H scattering vs. short distance 
cut-off Rmin- The laboratory kinetic energy is 3 eV. 


where the parabolic quantum number rii runs from 0 to 
n — \A\ — 1. Figure shows Vf,f{.{R) as a function of 
R for n = 3 and J = 4. Three of the nine potentials 
are attractive for small R and have a R~^ singularity in 
i? = 0. The corresponding phase shifts are ill-defined in 
these cases. When the correct angular coupling is used 
and the energy shifts are included, the potential curves 
are modified, but the problem with the ill-defined phase 
shifts in the dipole approximation remains. In the present 
model, the problem is cured by inserting an infinitely 
hard sphere with radius i?niin around the target nucleus. 
For low angular momentum, this introduces a dependence 
of the scattering matrix on the cut-off parameter Rmin- 
This should be considered as an uncertainty of the model 
related to the approximate treatment of the short dis¬ 
tance behavior. An example of the Rmin dependence in 
the calculation of the {Tr~p)n =3 + H scattering is shown 
in fig. for the /-average {I > 0) Stark, transport, and 
maximum absorption cross sections. Both the Stark and 
the absorption cross sections are rather insensitive to the 
value of i?min, whereas the transport cross section, which 
is more sensitive to the low partial waves, varies mod¬ 
erately (1.0 — 1.8 Ug). However, this dependence of the 
transport cross section is less significant when the energy 
distribution in the atomic cascade is taken into account. 


The statistically weighted differential cross sections for 
the I > 0 sector for (7r“p)„=3 J- H scattering shown in 
fig. 13 have the similar shape as those of the corresponding 
process for n = 5 in muonic hydrogen (fig. ^). Both the 
semiclassical and the fixed field model work well at kinetic 
energy T > 1 eV. 

The J dependence of the average Stark cross section 
(only the states I > 0 are included) and the maximum ab¬ 
sorption cross sections (the sum of the absorption during 
collision and 3p, 3d ^ 3s) at kinetic energy 10 eV is shown 
in fig. E- Like in the case of muonic hydrogen, the contri¬ 
bution to the Stark cross section can be divided into two 
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Fig. 13. Statistically weighted differential cross sections for 
the {'K~p)n ^3 +H scattering vs. CMS scattering angle 9 at the 
laboratory kinetic energy T = 10 eV. The solid line is the fully 
quantum mechanical result, the dashed line is the semiclassical 
model, and the dash-dotted line is hxed held model. 



Fig. 14. The average Stark and maximum absorption cross 
sections for the (7r“p)n=3 + H collisions vs. the total angular 
momentum J at the laboratory kinetic energy T = IQ eV. The 
solid lines are the fully quantum mechanical results, the dashed 
lines are the semiclassical model, and the dash-dotted lines are 
the hxed held model. The light solid line is the unitarity limit 
and the light dashed line corresponds to the assumption of 
statistical mixing. 


parts: a small J region, where the mixing is strong and 
the distribution over the final states is approximately sta¬ 
tistical, and a large J region, where the transitions with 
\Al\ = 1 dominate. The semiclassical model is in a good 
agreement with the quantum mechanical model, whereas 
the agreement with the fixed field model is fair. All three 
models are in a fair agreement for the maximum absorp¬ 
tion cross sections. 


The relative role of nuclear absorption during the col¬ 
lisions and between the collisions is illustrated in fig. |l^ 
for pionic hydrogen with n = 3. The absorption cross sec¬ 
tion (ie. the cross section for nuclear absorption during 
the collision, eq.(p^) is shown in comparison with the 
maximum absorption cross section, eq.(|3^. For energies 
larger than ^ 2 eV the tt~p atom is more likely to leave 
the collision zone in the 3s state than undergo nuclear 
absorption. At high density, many of the (7r~p)3s atoms 
with high kinetic energy will leave the 3s state before the 
nuclear reaction can take place. The results of the semi¬ 
classical and the fixed field model are in a good agreement 
with the quantum mechanical results for kinetic energies 
larger than 10 eV. Below 10 eV the semiclassical descrip¬ 
tion breaks down and significantly underestimates the ab¬ 
sorption cross sections. 

Figure shows the rates for different processes for the 
(Tr~p)n =3 state in hydrogen gas at a pressure of 15 bar. 
For energies larger than 0.05 eV Stark mixing is the fastest 
process while the 3p —> Is radiative transition dominates 
for low energies. The effective absorption rate is smaller 
than the Stark mixing rate due to the low statistical weight 
of the 3s state and the 3p — 3s energy difference (only im¬ 
portant for low energie^ The rates for 3^2 adiabatic 
Coulomb deexcitation |^] and external Auger effect 0] 
are also shown for comparison; these collisional deexcita¬ 
tion mechanism are obviously suppressed by the absorp- 
tion^. The deceleration rate exceeds the absorption rate 
below 2 eV, but for higher energies the deceleration is sup¬ 
pressed by the absorption, and for energies above 20 eV 
becomes insignificant. 


4.3 Kaonic hydrogen 

The nuclear interaction effects in the scattering of the 
-K~p atoms in excited states are even more important than 
in the case of 7r~p. The central values of the KEK result 
for the Is strong interaction shift and width of kaonic hy¬ 
drogen 1^ were used in our calculations: 

= 327±63(stat.) ± ll(syst.) eV , (80) 

Fi, = 407 ± 208(stat.) ± 100(syst.) eV . (81) 

The cascade in kaonic hydrogen differs from that of pi¬ 
onic hydrogen in the initial condition, it begins with a 
higher n level: ~ ^/pK-p ~ 25 as compared to 

'niniti'^~p) ~ y/I^TT-p ~ 15- The much larger width of the 
ns states makes the absorption during the collisions much 
more probable than in the 7r“p atom. Another important 
difference between the K~p and 7r“p cases is that in K~p 
the ns energy shift is repulsive, therefore the nl ns tran¬ 
sitions are not allowed below the corresponding threshold. 

^ There is experimental evidence for the 3^2 Coulomb¬ 
like deexcitation process in the neutron time-of-fiight spectra 
in both liquid and gaseous (40 bar) hydrogen which can 
be related to the formation of excited molecular states, see 
and references therein. 
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Fig. 15. The energy dependence of the absorption cross sec- 
tioiis, eq.(p^, and the maximum absorption cross sections, 
eq.(|3q), for pionic hydrogen with n = 3. The solid lines are 
the fully quantum mechanical results, the dashed lines are the 
semiclassical model, and the dash-dotted lines are the hxed 
field model. 



Fig. 16. The energy dependence of the /-average Stark 
(solid line), deceleration (dashed line), and effective absorp¬ 
tion (dash-dotted line) rates for (7r“p)n=3 at 15 bar. The light 
dotted lines are the radiative 3p —> Is and 3ci —> 2p rates. The 
light short-dashed line is the 3^2 Auger deexcitation rate Q 
and the light dash-dotted line is the 3^2 adiabatic Coulomb 
deexcitation rate 0 - 

The strong interaction width of the 2p stateQ is poorly 
known from the KN scattering data, but its effect on the 
collisional rates is negligible. 

Figure shows an example of the energy dependence 
of absorption cross sections for the states with n = 5. 
The quantum mechanical results are shown only above 
the 5s threshold as our numerical algorithm used in this 
particular case is not reliable in the presence of closed 

® It is important for the cascade calculations and strongly 
influences the yield of K X-ray lines, see |33[|34|]. 



Fig. 17. Absorption cross sections for K~p in the 51 states 
with Z = 1 — 4 scattering from hydrogen vs. laboratory kinetic 
energy. The fully quantum mechanical results are shown with 
solid lines, the semiclassical with dashed lines, and those of the 
fixed field model with dash-dotted lines. 

channels. Overall, there is a fair agreement between the 
three models above the threshold. The fixed field model 
does, however, result in absorption cross sections that are 
somewhat smaller for / > 1 than those of the semiclassi¬ 
cal model. This can be explained using hg. the correct 
angular coupling allows a larger fraction of the substates 
with / > 0 to be mixed with the s state in a single collision 
and thereby undergo nuclear absorption. Figure ^ shows 
the /-average Stark, absorption, and deceleration rates for 
n = 5 calculated in the semiclassical model in compar¬ 
ison with the deexcitation rates for a typical gas target 
at 10 bar. The {K~p)n =5 atoms with high kinetic ener¬ 
gies are strongly absorbed while the radiative and Auger 
deexcitations dominate in the low energy range. Like in 
the case of pionic hydrogen the situation is complicated 
by the deceleration due to elastic collisions. The results of 
cascade calculations based on the presented cross sections 
were discussed in [ p4[ . 

4.4 Antiprotonic hydrogen 

The case of antiprotonic hydrogen is similar to that of 
K~p: the ns nuclear widths are large and the ns nuclear 
shifts are repulsive. The following values |^] for the spin- 
averaged shift and width were used in the present calcu¬ 
lations: 

= 721 ± 14 eV , (82) 

As = 1097±42 eV . (83) 

The hadronic width of the 2p state in antiprotonic hydro¬ 
gen || 

= 32.5 ± 2.1 meV (84) 

is much larger than the radiative one and absorption dur¬ 
ing the cascade from the p states is very important. The 
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Fig. 18. The energy dependence of the /-average Stark, ab¬ 
sorption and deceleration rates for [K~p )„^5 at 10 bar. The 
results are calculated in the semiclassical model. The /-average 
Auger rates and the radiative rates are shown with light 
lines. 


Fig. 19. Absorption cross sections for pp in the 8/ states with 
/ = 1—7 scattering from hydrogen vs. laboratory kinetic energy. 
The results of the semiclassical model are shown with dashed 
lines, and those of the hxed held model with dash-dotted lines. 


(85) 


Figure |T^ shows the calculated cross sections for ab¬ 
sorption from the s state at n = 8. The energy depen¬ 
dence and the I dependence of these cross sections are 
very similar to the kaonic hydrogen case shown in fig. 0- 
The semiclassical and fixed field models are in a fairly 
good agreement. As was observed in the K~p case, the 
fixed field model underestimates the absorption cross sec¬ 
tions for I > 1 due to the approximate treatment of the 
angular coupling. The rates for the /-average collisional 
processes for {pp)n=8 at a pressure of 1 bar are shown 
in fig. ^ in comparison with the deexcitation rates: here 
the absorption just begins, it gets more important at the 
lower n, like in the {K~p)n^^ case shown in fig. 0 , and 
then eventually terminates the cascade. As the pp annihi¬ 
lation rate in the 8p state is 10^^ s“^, the p state is almost 
completely depleted between the collisions at 1 bar. The 
absorption between collisions from the 8p state is compa¬ 
rable in strength to the absorption from the 8s state dur¬ 
ing collisions for energies below 1 eV and about two times 
stronger for higher energies. The results of detailed cas¬ 
cade calculations for antiprotonic hydrogen will be pub¬ 
lished elsewhere. 


5 Conclusion 

The total and differential cross sections of Stark mixing 
and elastic scattering have been calculated for the p~p, 
7r“p, K~p, and pp atoms for the principal quantum num¬ 
bers and the kinetic energies needed in detailed cascade 
calculations. For hadronic atoms, the collisional absorp¬ 
tion cross sections have been calculated as well. For the 



Fig. 20. The energy dependence of the /-average Stark (solid 
line), effective absorption (dash-dotted lines) and deceleration 
(dashed line) rates for pp for n = 8 at 1 bar. The rates are calcu¬ 
lated in the semiclassical model. The /-average Auger rates 
and the radiative rates are shown with the light lines. 


low states n = 2 — 5, the calculations have been done in 
a fully quantum mechanical framework using the close¬ 
coupling method. For the first time, the effects of nuclear 
shifts and width of the ns states have been taken into ac¬ 
count straightforwardly in the quantum mechanical scat¬ 
tering problem. For the intermediate states n = 5 — 10, the 
proposed semiclassical model provides an efficient compu¬ 
tational method. By treating one degree of freedom (the 
distance between x~p and H) classically, one can reduce 
the original system of the coupled second order equations 
to a system of first order equations while maintaining the 
correct angular coupling between the x~p internal angu¬ 
lar momentum and the relative orbital angular momentum 
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of x~p + H system. The semiclassical approximation has 
been found to agree fairly well with the fully quantum 
mechanical calculations, provided the collisional energy is 
not too low so that the number of essential partial waves is 
large. As the threshold behavior is not treated correctly in 
the semiclassical approximation, transitions to and from 
the ns states cannot be calculated reliably in the near¬ 
threshold region. 

Using the above described methods we were able to 
assess the range of validity of the fixed field model, which 
was commonly used in many cascade studies. In addition 
to total Stark mixing and absorption cross sections cal¬ 
culated with this model in the literature 00 ] , we have 
calculated the differential cross sections for comparison 
with the more accurate methods. When compared to the 
semiclassical model, the fixed field model usually under¬ 
estimates the absorption cross sections due to the lack 
of the rotational coupling among the nl sublevels. When 
nuclear absorption during collisions is negligible, the fixed 
field model provides, on average, a fair description in com¬ 
parison with the semiclassical approximation for kinetic 
energies larger than a few eV. 

The results of this paper have been used in detailed 
kinetics calculations of atomic cascade in fj,~p and 'K~p 
reported in [^. The detailed description of the results, 
together with those in K~p and pp, will be published in 
separate papers 

A few problems remain, which are beyond the scope of 
this paper. First, we have considered collisions with atomic 
hydrogen. One can expect that the molecular structure of 
the target becomes important for large n states when the 
characteristic size of the exotic atom cannot be treated 
as a small parameter in comparison with the conventional 
atomic scale. We shall address this problem in a sepa¬ 
rate paper [Q . Molecular effects @ are also expected to 
be important at low collisional energy when only a small 
number of molecular ro-vibrational states can be excited. 
This kinematical region partly overlaps with the region of 
small energies where only a few partial waves are impor¬ 
tant and the dipole approximation is not justified. To deal 
with these problems a genuine many-body framework is 
needed. 
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A The rotation functions {3,j) 

This paper uses the conventions of Condon and Short- 
ley for the spherical harmonics and those of ref. |pM for 
the rotation functions Let I? = (0,0) and 

C' = (0', 0') be a direction expressed the original and the 
rotated coordinates, respectively. The effect of the rotation 
on a spherical harmonic is given by the rotation functions 




YjA{n') = E Y>ma ■ (Al) 


M 


The rotation functions Pj'j) are eigenfunctions of 

the square of the total angular momentum and its pro¬ 
jections Jz and Jz' along the 2 ;-axis and the z'-axis: 


= J{J+ , 

JzDiiA{a, /?, l) = f3, j) , 

Jz'DiiAia, 13, 7) = -ADij^ia, 13, 7) . 

The raising and lowering operators are defined by 

— Jx ^ '^Jy 1 

(/j_ = Jx' i iJy' 

and have the properties 

J±DMA{a,(3,'y) = -X±{J, A)DiiA±i{a, (3,-i) 


where 


\±{J,M) = yj(jTi)^wu(M^T) 


(A2) 


(A3) 


(A4) 


(A5) 


In this paper only two of the Euler angles, a and (3, are 
used. Therefore, to simplify notation the rotation func¬ 
tions can be written ( 12= (0, 0)): 


DitAm = DitA{^,e,Q) 


(A6) 


The following integrals represent the normalization and 
orthogonality conditions for the rotation functions and the 
Wigner-Eckart theorem for the irreducible representations 
of the rotation group 


/dI2 D{^Ai^)Y>M'A'{^) — 


5jj'5aa'3mm' 


Jtt 


2J-k 1 ’ 

Jtt 


(A7) 


2J+1 


{LMLlm\JM){JM'\LM'^lm') (A8) 


where {LM]^lm\JM) denotes a Clebsch-Gordan coefficient. 
For A = 0, the rotation functions are related to the spher¬ 
ical harmonics by 


r)J 

^MO 


( 12 ) = 


47r 


2J+1 


-V 


jm(^) 


(A9) 


The coefficients Ual use d in the ba sis t ransformation 
can be found using eqs. (AS) and (A9) 


^AL — 


2J+1 


47r 


J dHdr (^Dil{Q)xniA (rO) T/f" (12, a;)r„i (r) 
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= (^{LlMLm\JM) 

Mliti 

X J dn DiM{[2)YLML{^) J drx*niA{r')Xnifn{r)'j 

= Y [{LlMLm\JM) 

MLm 

X J dn 

= ^J^^{L10A\JA) (AlO) 

for any M with \M\ < J. 

B Variable phase approach to multichannel 
scattering 

We use a version of the variable phase method (see ref. 
and references therein) to compute the scattering matrix. 
In this approach the problem of solving the second order 
linear and homogeneous Schrodinger equation is trans¬ 
formed into a problem of solving a nonlinear first order 
equation for the scattering matrix. 

In this appendix the notation is as follows: ^ is a col¬ 
umn vector containing the radial wave functions, K is 
the diagonal matrix with the channel momenta Kmn = 
dmnkm- The angular momentum quantum numbers for the 
different channels are diagonal elements of the matrix L, 
i.e. Lmn = dmrdm- Then the radial Schrodinger equation 
is given by 


Foll owing the same procedure as in the derivation of eq. 
(B2) one obtains 


S' = - H2)Weff{HiS - H 2 ) (B4) 

where Hi{R) = and H 2 {R) = . 

The matrix S{R) is related to the scattering matrix by 


S = . (B5) 

For short range potentials eq. appears to be less suit¬ 
able than eq. ( p^ ) because the solution must be obtained 
for large values of R due to the long range behavior of 
centrifugal repulsion. This is, however, not necessary be¬ 
cause it is always possible to convert S{R) to S{R) and 
vice versa with algebraic methods. 

The connection between S{R) and S{R) can be es¬ 
tablished through the identity of the wave functions and 
their derivatives in R. Let S' be a square matrix with lin¬ 
ear independent solutions ^ as colujnns. S and S' can be 
expressed both in terms of S and B 


S' = {HiS - H2)N = {IRS - H2)N 
S'^ {H'^S - H'^)N = {H'-^S - H'^)N (B6) 


where N and N are square matrices. From the relations 
(B6), one finds the following expression for S 


s = h: 


HiS - H2)NN-^ + H 2 


(B7) 


with 


NN-^ = [h'^{HiS - H2) - Hi{H'^S - 

X (iJiiJ' -H'^H2) . (B8) 


where W{R) is the reduced potential matrix. 


Let and be Riccati-Hankel functions as de¬ 
fined in ref. |^| and Hi and H 2 diagonal matrices with el¬ 
ements R')/ "x/and H 2 mn{R{ — 

r)/Vk^. 

The scattering matrix S{Rq) obtained from eq. ( pl| ) 
with W truncated at Rg {i.e. with the substitution W 
W{R)6{Rq — R)) is a function of Rg and satisfies the equa¬ 
tion 


,(2) 


= -{SHi - H2 )W{HiS - H 2 ) 


(B2) 


The scattering matrix of the full problem is given by S' = 
S(oo). 

In eq. (p32| ) the dependence on the angular momentum 
is contained in the Riccati-Hankel functions. In numerical 
calculations, it can be more convenient to combine the 
potential and the angular momentum term in one effective 
potential 


+ . (B3) 


We compute the scattering matrix for the x~p + H ^ 
x~p + H process by solving eq. ( |b^ in the coupled basis 
with the boundary condition that S(Rmin) is a diagonal 
matrix with elements 


= Sr, 


1 F ^ tan(/i:^Rjnin) 
1 i tan(A:^Rjnin) 


(B9) 


which is the scattering matrix for S-wave scattering from 
an infinitely hard sphere with radius Rmin- The effect 
of nuclear absorption from the ns states is included by 
adding the imaginary part (the real part is already taken 
into account in the momentum matrix K) of the ns energy 
shift to the potential 


IFy(R) = 2p{{n]L'l'JM\V{R)\n-LlJM) 

Y5ig5vg{-irrrsl2)e{Rg-R)) . (BIO) 

The nuclear absorption is turned off for distances larger 
than Rg\ we use Rg = 5ag in this paper. This allows the 
mixing nl <-*■ ns to take place during the collision together 
with the absorption effects while the ns states remain well 
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defined asymptotic states. The nuclear reactions that oc¬ 
cur after the collision are taken into account by a cascade 
model. 

It is necessary to set the boundary condition away 
from R = 0 because the potential taken in the dipole 
approximation has an R~^ singularity at i? = 0 which 
makes the Schrbdinger equation ill-defined. The exact po¬ 
tential, however, becomes finite for small R. By using 
the short distance cut-off i?min we obtain a well-defined 
Schrodinger equation, but the calculated cross sections 
become dependent on this regularization parameter. By 
varying i?niin around the value where the exact potential 
becomes repulsive one can get an estimate of the uncer¬ 
tainty of the cross sections due to the treatment of the 
short distance behavior. Throughout this paper we use 
the value i?min = 0.05ao unless otherwise stated. 
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